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Abstract 

In order to determine polarized parton densities in nucleon we have 
made fits using all experimental data on spin asymmetries measured 
in the deep inelastic scattering on different nucleon targets. We have 
used in our analysis next to leading order QCD corrections in MS 
renormalization scheme. The functional forms of polarized densities 
are based on the fit to unpolarized deep inelastic data made by the 
MRST group (MRST 99). We have concentrated on different models 
for gluon spin distributions. In the best fit we get for gluon polariza- 
tion (at Q 2 = lGeV 2 ) AG = -1.0 ± 0.6. The obtained result does 
not depend very much on various model assumptions. 
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The experimental data on deep inelastic scattering of leptons on nucleons 
enable us to determine parton densities in nucleons. For polarized particles 
such experiments were made in SLAC 0, CERN § and DESY §. The 
results were analyzed by many groups and several spin distributions and 
different gluon polarization were presented [§, |[ ||, 0. In this paper we 
do not consider any new experimental data. All existing data on polarised 
deep inelastic scattering were already included in our latest fit where we 
have determined polarized parton densities (performing next to leading order 
QCD analysis in MS renormalization scheme). However, since our method 
of determination of spin distributions depends very strongly on the parton 
distributions for unpolarized case we have performed a new analysis using 
unpolarized distributions obtained by Martin, Roberts, Stirling and Thorne 
(MRST) @. 

The unpolarized densities are the sums of + and — components (densities 
of quarks and gluons with helicity along or opposite to the helicity of the 
parent nucleon) whereas the polarized ones are given by the differences. The 
asymptotic behaviour of our polarized parton distributions is determined (up 
to the condition that the corresponding spin densities are integrable) by the 
fit to unpolarized data. 

The determination of unpolarized parton densities performed by MRST 
[0] gives three solutions for gluon distribution. In this paper we want to follow 
the method used in our previous papers || [7| and to check how different 
functional form of gluon distributions influences our fits. We will use data 
for spin asymmetries at given x and different Q 2 . Altogether we will use in 
our analysis 431 experimental points. In our formulae for polarized densities 
we have included 12 parameters. 

Experiments on unpolarized targets provide information on the spin av- 
eraged quark and gluon densities q(x, Q 2 ) and G(x, Q 2 ) inside the nucleon. 
In ref.|| these distributions (at Q 2 = 1 GeV 2 ) are parametrized as follows: 

qi (x) = AiX X *(l - x) m {l + + fax), (1) 

where one assumes that sea antiquarks and quarks of the same flavour have 
identical distributions. The similar form one has for gluon density: 

G(x) = B g x Xg (1 - xf G (l + e GV ^ + nax). (2) 

The values of constants A» and rji, used in eqs. (1) and (2), are given in Table 
1 for all three solutions which are called lower (L), central (C) and higher 
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(H) gluon solution in ref.H (the differences are visible for big x values). Our 
polarized densities for valence quarks (at Q 2 = 1 GeV 2 ) are: 

Aqi(x) = x A *(l - x) Vi (ai + fayfx + 7i x), (3) 

whereas for sea quarks and antiquarks one has: 

Aqj(x) = x x '^{l-x) ns {a j +l3 j y/x)+CiX Xs {l-x) ns+2 as{l+es^/x+iisx), (4) 

where q = +1(— 1) for up (down) flavour and zero for strange quarks. With 
such parametrization of quark spin densities we get from the fit approxi- 
mately as ~ and a-u ~ otq ~ 2«s > Ar ^ Hence, our a posteriori 
assumptions are: 

a 5 = 0, 

olu = ocj = 2a-s, (5) 
At = Aj- 

Such choice is consistent with the SU(2) symmetry in the sea. 
For polarized gluon distribution we assume: 

AG{x)=x XG {l-xf G {a G + (3 G ^ + lG x). (6) 

The parameters a^ fa and 7j are determined from our fits and presented in 
Table 2. 

Table 1. The values of constants which appear in eqs. (1-4,6) for three 
solutions of unpolarized gluon density. This figures determine the behaviour 
of distributions at x = and x — 1 . 



Constant 


Lower 


Central 


Higher 




gluon solution 


gluon solution 


gluon solution 




-0.527 


-0.583 


-0.560 


Vu v 


3.44 


3.43 


3.45 


^d v 


-0.748 


-0.730 


-0.740 


Vd v 


3.86 


3.86 


3.95 


K 


-1.240 


-1.282 


-1.273 


Vs 


7.63 


7.65 


8.30 


\s 


-0.356 


0.183 


0.157 


A G 


0.131 


-0.031 


-0.005 


Vg 


5.67 


6.88 


7.45 
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Table 2. The parameters of three fits calculated at Q 2 = 1 GeV 2 . 



Parameter 


Lower 


Central 


Higher 




gluon solution 


gluon solution 


gluon solution 




0.110 


0.499 


0.279 




-3.50 


-4.71 


-4.22 




13.3 


14.0 


13.9 


«*, 


-0.055 


0.028 


-0.059 


/?*, 


-1.66 


-1.70 


-1.67 


1d v 


0.007 


-0.077 


-0.007 


Ctu 


-0.033 


-0.133 


-0.082 


At 


1.54 


1.37 


1.58 




-0.402 


-0.056 


-0.226 




-18.3 


-24.1 


-20.3 


Pg 


52.5 


117.2 


86.7 


lG 


-48.8 


-173.5 


-127.9 



The total quark distributions are given by: 

Au = Au v + 2A% 

Ad = Ad v + 2Ad, (7) 
As = 2As, 

whereas one has for axial charges: 

AS = Au + Ad + As, 
a 8 = Au + Ad-2As, (8) 
a 3 = Qa = Au — Ad. 

In Table 3 one can find above quantities integrated over variable x and 
calculated at Q 2 = 1 GeV 2 . The axial charge a% is not fixed (usually one 
fixes it in an analysis) and comes out close to its experimental value. The 
quantity a 8 is in some sense fixed i.e., it is included in a fit as an additional 
experimental point. We put a s = 0.58 ±0.1, the value taken (with enchanced 
to 3<r error) from the data on hyperon (3 decays. 

The typical errors of our determination are: AG = —1.02 ± 0.61, AE = 
0.04 ± 0.18, where we quote the figures for higher gluon solution. We use all 
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three MRST solutions keeping the quantity as fixed. As it is seen from Table 
3 the best x 2 we g°t for higher gluon solution. 

Table 3. The integrated quantities from our three fits calculated at Q 2 = 
lGeV 2 . 



Quantity 


Lower 


Central 


Higher 




gluon solution 


gluon solution 


gluon solution 


AS 


0.19 


0.02 


0.04 


a 8 


0.58 


0.59 


0.51 


9a 


1.25 


1.29 


1.31 


Au 


0.79 


0.75 


0.75 


Ad 


-0.46 


-0.54 


-0.55 


As 


-0.13 


-0.19 


-0.16 


AG 


-0.77 


-1.15 


-1.02 


x 2 


363.91 


363.16 


362.28 



For higher gluon solution we have tried to make several fits, where we 
investigate other assumptions which slightly change gluon distribution. In 
one of the fits we put chq h =0 at Q 2 = 1 GeV 2 (it means, that gluon density is 
less divergent at x = 0). In this case we get a fit which is of undistinguishable 
quality from the best fit (in both cases \ 2 P er degree of freedom is 0.86) and 
one gets AG = - 0.63. However, when we put ac H = (3g h = 1g h = (i.e. 
AG = at Q 2 = 1 GeV 2 ) we get x 2= 366.0 (per degree of freedom we have 
0.87) which gives not much worse fit. We have also performed other fits, 
one in which a§ is treated as free parameter and in another we allow gluon 
distributions G + (x) = (G(x) + AG(x))/2 and G~{x) = (G(x) - AG(x))/2 
not to be positively defined. In both cases we have AG = —1.0 and \ 2 is 
nearly as good as in the best fit. In all considered models we get negative 
polarization of gluons (the negative value of gluon polarization one can find 
also in ||). The values of gluon polarization (from our best fit) at other Q 2 
values are: AG(Q 2 = 5 GeV 2 ) = -2.0 and AG(Q 2 = 100 GeV 2 ) = -3.9. 
The variations of our model does not change significantly other quantities 
(see Table 3). 

In figures 1 and 2 we present the shapes of gluon spin distribution obtained 
in our fits at Q 2 = 1 GeV 2 . One can see that the curves does not differ very 
much. 
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Figure 1: The comparison of polarized gluon densities xAG(x) versus x ob- 
tained from the fits corresponding to lower gluon (dotted line) central gluon 
(dashed line) and higher gluon (solid line) solutions from ref.MI. 
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Figure 2: The comparison of polarized gluon densities xAG(x) versus x ob- 
tained from the fits corresponding to the higher gluon (solid line) and the 
same model with cig h = (dashed line). 
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We have made fits to precise data on spin asymmetries on proton, neutron 
and deuteron targets. Our model for polarized parton distributions is based 
on MRST 1999 fit || to unpolarized data. Gluon polarization comes out 
negative and relatively high. This conclusion does not depend very much on 
different assumptions (within our famework) about gluon spin distribution. 

PACS numbers: 12.38.-t, 13.60.Hb, 13.88. +e, 14.20.Dh 
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